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MEASURED AND CALCULATED FAST NEUTRON SPECTRA IN A DEPLETED 
URANIUM AND LITHIUM HYDRIDE SHIELDED REACTOR 
by Gerald P. Lahti and Robert A. Mueller 
Lewis Research Center 


SUMMARY 

Measurements of MeV neutron spectra were made at the surface of a cylindrical 
shield arrangement surrounding a small (25. 4 -cm diam), cylindrical uranyl fluoride 
(UOgFg) - water solution reactor. Four shield configurations were considered: These 
were assembled consecutively with cylindrical shells of 5 -centimeter -thick depleted ura- 
nium, 13 -centimeter -thick lithium hydride, 5 -centimeter -thick depleted uranium, and 
3 -centimeter -thick depleted uranium. Measurements were made using a NE-218 scin- 
tillation spectrometer. Proton recoil pulse height distributions were differentiated to 
obtain neutron spectra. Calculations of leakage flux were made with the two-dimensional 
discrete ordinates code DOT using ENDF version 3 cross sections. Comparisons of ex- 
perimental and calculated results on an absolute basis show some deviations, with the 
deviations growing with increasing shield thickness. This discrepancy may be due to 
cross-section errors. However, excellent agreement between measured and calculated 
spectral shapes was observed. 


INTRODUCTION 

The design of a minimum weight shield for a nuclear power reactor requires con- 
fidence in calculations and cross-section data. To check neutron transport calculations, 
a number of transmission experiments were conducted using both small neutron sources 
and the NASA Zero -Power Reactor Facility. This series of experiments first examined 
geometrically clean spherical configurations of lead and other materials (refs. 1 and 2) 
and then cylindrical configurations of lead and water (ref. 3). 

The present report describes calculations and measurements made for the cylin- 
drical reactor surrounded by alternate layers of depleted uranium and lithium hydride. 
Calculations were made using the two-dimensional discrete ordinates transport code 



DOT (refs. 4 and 5). Neutron cross sections were from ENDF/B (version 3) files; the 
basic data were collapsed to a 4 3 -group multigroup library using SUPERTOG-1H (ref. 6). 
Measurements were made using a 1. 22 -centimeter inside diameter by 1. 27-centimeter 
long NE-218 scintillator. The measurements of proton recoil spectra were differentiated 
using the TUNS code (ref. 7) to obtain neutron spectra. 


EXPERIMENTAL SETUP 
Facility 

The test facility used in these experiments is the NASA Zero Power Reactor Facil- 
ity. The uranyl fluoride - water solution reactor in this facility is 25.4 centimeters in 
diameter, and, depending on the fuel concentration and reflector material, can be oper- 
ated at heights up to 76 centimeters. The reactor is licensed to operate at a maximum 
of 100 watts thermal (total neutron source strength of about 8x10* neutrons/sec). The 
reactor vessel is 0. 278 centimeter thick and is made of AISI 316 stainless steel. The 
reactor is surrounded by a 244 -centimeter -diameter by 104 -centimeter -deep cylindrical 
tank; shielding materials may be arranged about the reactor in this tank. 


Present Experimental Setup 

The present experimental setup is illustrated schematically in figure 1. Exact di- 
mensions and material descriptions are given in table I. The reactor is surrounded by 
91 -centimeter -high cylindrical shells of water (5 cm), depleted uranium (5 cm thick), 
lithium hydride (13 cm thick) followed by two more of depleted uranium (5 cm and 3 cm 
thick). Midplane measurements of the neutron transmission were made for these four 
configurations: 

(1) Core plus water plus the first 5-centimeter -thick depleted uranium layer 

(2) Configuration 1 plus a 13 -centimeter -thick lithium hydride layer 

(3) Configuration 2 plus a 5 -centimeter -thick depleted uranium layer 

(4) Configuration 3 plus a 3 -centimeter -thick depleted uranium layer. 

The thicknesses of these layers correspond roughly to those considered for space power 
reactor shields. The midplane measurements were made with the detector at the outer 
surface of each configuration. 
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Shield Material Fabrication 


The depleted uranium shells were rolled from plate and canned with 0. 19 centimeter 
(75 mil thick) stainless steel. The lithium hydride shells were cast in segments using 
techniques developed in connection with the Atomic Energy Commission's SNAP reactor 
program. The lithium hydride was internally reinforced using a stainless -steel honey- 
comb matrix. About a 94 -percent theoretical density of lithium hydride was achieved in 
those castings. Details of this lithium hydride casting may be found in reference 8. 


MEASUREMENTS 

Spectrometer 

Measurements of the neutron spectra emerging from the shield assembly were made 
at the horizontal midplane using a 1. 22 -inside -diameter by 1. 27 -centimeter -long NE-218 
scintillator. The scintillator was optically coupled to an RCA 8575 photomultiplier tube 
using a 1. 3 -centimeter light pipe. A modified Owens type of pulse shape circuit (as de- 
scribed in ref. 9) was used. This circuit is used with a two-parameter pulse-height an- 
alyzer to separate proton recoil (neutron -induced scattering events) and electron recoil 
(gamma -ray -induced scattering events) occurring within the scintillator. A spectrometer 
block diagram is illustrated in figure 2. The proton recoil region of the two -parameter 
analyzer data was summed using the PREJUD code (ref. 10); the output data from 
PREJUD was then differentiated using the TUNS code (ref. 7) to obtain neutron spectra. 
Reference 7 also contains more details on this experimental technique. 


Power Calibration 

To compare calculated results with the experimental measurements on an absolute 
basis, both were normalized to correspond to a reactor thermal power level of 1 watt. 
The reactor power was calibrated by mapping the thermal flux within the reactor with 
dysprosium foils, which had previously been calibrated against indium foils. The total 
fission rate was then calculated from the thermal neutron flux level determined by these 
foil measurements. This was then correlated with integral counting data taken over the 
period of measurement using a helium -3 proportional counter. The helium -3 propor- 
tional counter data was then used to determine reactor power levels for all other 
measurements. 
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Experimental Error 


Systematic errors . - The uncertainty in the absolute power level of the reactor is 
estimated to be ±15 percent. Because of the nature of the method of determining power 
level, this relative error will be approximately constant throughout the present series of 
experiments and is one possible source of systematic error. The relative power level 
uncertainty is less than 5 percent. 

Another source of systematic error is introduced through the detector efficiencies 
used in the spectrum unfolding code. The efficiency data used was that for a slightly dif- 
ferent detector and was measured without a light pipe between the detector and the photo- 
multiplier tube. One test measurement, however, showed little difference in efficiencies 
between the detector systems; thus, only a small error is anticipated from this 
uncertainty. 

Random errors . - Counting statistics introduce a random error to the measure- 
ments. These are represented by the error bars appearing in the plots of experimentally 
measured neutron spectra and are an inverse measure of the count rate observed. 


CALCULATIONS 

The leakage flux per unit energy in the four experimental configurations described 
earlier was calculated using the two-dimensional discrete ordinates code DOT. DOT was 
run in a reactivity -cal culation mode to account for neutron production in the outer ura- 
nium layers. A symmetric S-8 (48 angle) quadrature and P-3 scattering approximations 
were used. The experiment was mocked up for calculation as described in table I. 

Table I also lists the number of radial mesh intervals considered. Midplane symmetry 
was assumed; 27 axial intervals were used in the calculation. Table n lists boundaries 
for the 43 -group split used in the calculations. The top 30 groups correspond to a leth- 
argy split of Au = 0. 1. The resonance and thermal groups were necessary in the cal- 

ooc ooo 

culation because of the neutron production in both 0 and U in the depleted uranium 
shield layers. 

Multigroup cross sections used were assembled from ENDF/B (version 3) data using 
the SUPERTOG-DI code. Data for several elements were taken from GAM-II files. The 
sources of cross-section data used is listed in table m. Macroscopic cross sections 
were assembled using atom densities listed in table IV. 
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COMPARISON OF MEASURED AND CALCULATED SPECTRA 


The neutron leakage flux from each of the four experimental configurations is shown 
in figure 3. Flux per unit energy is shown for the interval of 0. 9 to 14 MeV. To provide 
a comparison on a uniform basis, both calculated and measured results were normalized 
to correspond to a power level of 1 watt (reactor fission rate of 3. ixlO 10 fissions/sec). 

Several stages of postprocessing were necessary to convert the DOT -calculated 
multigroup fluxes to those shown in figure 3. First, the surfaces angular fluxes were 
converted to a scalar flux at the surface using the equation 

<p= XI w i^i 

All 

angles 

where (p is the group scalar flux, w. is the weight (solid angle) associated with if /. , and 

T fh 1 1 

i//. is the i angular flux. Next, this surface multigroup scalar flux was converted to 
flux per unit energy E by dividing it by the energy group width, that is, 


<p{ E) = 

AE 

Finally, because the optics and the electronics of the instrumentation system introduce 
statistical variations, a Gaussian smearing of data is performed by the instrumentation. 
To provide a realistic comparison between calculation and experiment, the calculated re- 
sults must also be smeared by this Gaussian function. Symbolically, this is 

s+OO 

cp s (E)dE = <p(E , )G[R(E'), E’, E]dE* dE v 

where ^E'ME' is the actual flux of neutrons in dE* about E T ; G[R(E’),E',E]dE is the 
probability (a Gaussian) that a particle of energy E' will appear to the instrumentation 

as a particle in dE about E, where G is of the form A/R(E’) exp (-B { (E' - E)/ 

[R(E')] } )^ and A and B are constants such that ^°° G dE = 1, namely, A = Vln 2/77 
and B = In 2; R(E') is the resolution (more precisely, the half width at half maximum) 
of the instrumentation for particles of energy E'; and <p (E)dE is the "smoothed" neu- 
tron flux in dE about E, which should be comparable to what the instrumentation sees. 

The resolution R(E) of the present instrument is taken to be half that used for the 
5_ by 5-centimeter (2- by 2-in. ) NE-213 detector (ref. 11). Only the smeared calculated 
fluxes are shown in figure 3. 
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In all cases, the measured and calculated spectral shapes (fig. 3) show excellent 
agreement. Differences in absolute magnitude greater than the uncertainty in the power 
calibration appear in the results for configuration 4. There also appears to be an in- 
crease in the difference between measured and calculated results with increasing shield 
thickness. Some possible reasons for the differences include: 

(1) Cross sections in error. A 5 percent error in, say, the total cross section of 
depleted uranium would result in a 20 percent error in calculated transmission through 
the 13 centimeters of depleted uranium considered here. Because differences are 
greater when lithium hydride is present, lithium hydride cross sections are more suspect. 

(2) Cracks and voids in the cast lithium hydride (as described in ref. 8), which could 
permit higher transmission than calculated for homogeneous material. 

(3) Inaccurate transport calculations due to inadequate spatial and angular mesh and/ 
or scattering order. But this is presently discounted; an S angular quadrature was 
used in one calculation and showed only a few percent difference from the Sg 
calculations . 

The increasing discrepancy with added thickness can most reasonably be attributed 
to cross section errors. Because of the good agreement in spectral shape, the inelastic 
cross section is not suspect. This leaves the total and/or elastic scattering cross sec- 
tion in doubt. But a definitive determination of such a discrepancy cannot be determined 
from an integral experiment such as this. 

On a positive note, the calculated results are all within 50 percent of the measured 
results. Thus, an engineering design tool, the calculative method (and cross sections) 
presently used is accurate to within this value. 


CONCLUDING REMARKS 

Measurements and calculations of MeV neutron spectra emerging from a depleted 
uranium and lithium hydride shielded reactor have been made. In general, the agree- 
ment between experiment and calculation is good. Spectral shapes are in excellent 
agreement. Comparisons made on an absolute basis are within 50 percent for all cases 
considered. Differences in spectra of the order of 50 percent were observed in the con- 
figuration using the thickest shield, namely, 13 centimeters of depleted uranium and 13 
centimeters of lithium hydride. ' This difference is possibly attributable to inhomogene- 
ities (cracks and voids) in the cast lithium hydride or to erroneous cross sections. 
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Latest ENDF (version 3) cross sections for lithium 6 and 7, hydrogen, uranium 235 and 
238, and iron were used in the calculations. 


Lewis Research Center, 

National Aeronautics and Space Administration, 
Cleveland, Ohio, July 2, 1973, 

503-05. 
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TABLE n. - ENERGY GROUP SPLIT USED IN CALCULATIONS 
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TABLE HI. - SOURCES OF CROSS 


SECTION DATA USED IN 


PRESENT CALCULATION 


Element 

ENDF material number 

Hydrogen 

1148 

Lithium -6 

1115 

Lithium -7 

1116 

Oxygen 

1134 

Fluorine 

(GAM-II data used) 

Chromium 

(GAM -II data used) 

Iron 

1180 

Nickel 

(GAM -II data used) 

Uranium‘^ i) 

1157 

Uranium 2 ’*** 

1158 


TABLE IV. - ATOM DENSITIES OF MATERIALS 


USED IN PRESENT CALCULATION 


Material 

Element or 
isotope 

Density, 

atoms/cm'* 

Core 

Hydrogen 

Oxygen 

Fluorine 

Uranium ^ 

IT . 238 

Uranium 

0. 06639X10 24 
. 033419 
. 00022236 
. 00010362 
. 00000756 

Stainless steel 

Chromium 

Iron 

Nickel 

0. 01559X10 24 
. 060638 
.009745 

Lithium hydride 

Hydrogen 
Lithium -6 
Lithium -7 

0. 056837X10 24 
.004217 
. 052620 

Depleted uranium 

Uranium 22 ** 

TT . 238 

Uranium 

94 

0. 000109X10 
. 047200 





r Reactor vessel 316 stainless steel 


t 


t 



Figure 1. - Sketch of experimental setup. Region above horizontal midplane is shown. Dimensions 
indicated are approximate and are only for general orientation. The exact dimensions are given 
in table I. 


RCA8575 



voltage 

Figure Z - Spectrometer block diagram and photomultiplier tube base circuit 
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Energy, E, MeV 

(c) Shield configuration: 5 centimeters of depleted uranium plus 13 centimeters (d) Shield configuration: 5 centimeters of depleted uranium plus 13 centimeters 

of lithium hydride plus 5 centimeters of depleted uranium. of lithium hydride plus 5 centimeters of depleted uranium plus 3 centimeters 

of depleted uranium. 

Figure 3. -Concluded. 
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